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Many meson processes are related to the Ua( 1) axial anomaly, present in the Feynman 
graphs where fermion loops connect axial vertices with vector vertices. However, the cou¬ 
pling of pseudoscalar mesons to quarks does not have to be formulated via axial vertices. 
The pseudoscalar coupling is also possible, and this approach is especially natural on the 
level of the quark substructure of hadrons. In this paper we point out the advantages of 
calculating these processes using (instead of the anomalous graphs) the Feynman graphs 
where axial vertices are replaced by pseudoscalar vertices. We elaborate especially the 
case of the processes related to the Abelian axial anomaly of QED, but we speculate that 
it seems possible that effects of the non-Abelian axial anomaly of QCD can be accounted 
for in an analogous way. 
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1. Introduction 

Numerous pr ocess es in meson physics are related to the Adler-Bcll-Jackiw (ABJ) 
axial anomaly ^ 21 appearing in the fermion loops connecting certain number of axial 
(A) and vector (V) vertices. Concretely, in this paper we will deal with the pro¬ 
cesses related to the AW (“triangle”, Fig.Q and VAAA (“box”, Fig. [2j anomaly, 
exemplified by the famous 7 r° —> 77 and 7 —> 7r + 7r°7r - transitions. 

Suppose one wa nts t o describe such processes using QCD-related effective chiral 
meson Lagrangians — — without adding ad hoc interactions of mesons with external 


* Senior Associate of Abdus Salam ICTP 


1 








2 


D. Kekez, D. Klabucar and M. D. Scadron 



Fig. 1. The triangle graph and its crossed graph relevant for the interaction of the neutral pseu¬ 
doscalar meson of momentum P with two photons of momenta k and k '. The quark-photon coupling 
is in general given by dressed vector vertices T^( 91 , 92 )? which in the free limit reduce to eQq^. 



Fig. 2. One of the box diagrams for the process 7 —» 7 r'* _ 7 r 0 7 r — , studied by the pseudoscalar 
coupling method in, e.g.,Refs. 3, 4. There are six different contributing graphs, obtained from the 
above graph by the permutations of the vertices of the three different pions. The position of the u 
and d quark flavors on the internal lines, as well as Q u or Qd quark charges in the quark-photon 
vertex, varies from graph to graph, depending on the position of the quark-pion vertices. The 
physical pion fields are 7r± = ( 7 ri =p 27T2)/\/2 and 7 r° = 7 r 3 . The momenta flowing through the 
four sections of the quark loop are conveniently given by various combinations of the symbols 
/?, 7 = +> 0 , - in k a/37 = k + api + (3p 2 + 7 P 3 - 


gauge fields to reproduce empirical results. For example, one can add by hand 

AC = , ( 1 ) 

and this would reproduce the observed n° —> 77 width for the favorable value of 
the 7r°77 coupling g n77 . However, if one does not want to add such ad hoc terms 
in the effective meson Lagrangians, one must describe such “anomalous” processes 
through the term derived by Wess and Zumino ^. On the other hand, if one wants 
to utilize and explicitly take into account the fact that mesons are composed of 
quarks, another way of describing these processes is optimal in our opinion, and the 
main purpose of this paper is to stress and elucidate this. 

Axial vertices in the anomalous graphs such as the AVV and VAAA ones, couple 
the quarks with pseudoscalar mesons. Instead of anomalous graphs, another way 
to study the related amplitudes involving pseudoscalar mesons, is to calculate the 
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corresponding graphs where axial vertices (A) are replaced by pseudoscalar (P) 
ones. Thereby, for example, the 7r° —► 77 decay amplitude due to the AVV “triangle 
anomaly”, 

° 2 N r 


F m „= o(tt u -> 2 7 ) = 


( 2 ) 


is reproduced by the calculation of the PVV triangle graph. [Eq. © pertains to the 
chiral limit, where the pion mass m^ = 0 . Also, ~ 93 MeV is the pion decay 
constant, e is the proton charge, and N c = 3 is the number of quark colors.] 

The PVV triangle graph calculation of Eq. 0 can most simply be done essen¬ 
tially a, la Steinberger ™ that is, with a loop of “free” constituent quarks with the 
point pseudoscalar coupling (i.e., 375, where g = constant) to quasi-elementary pion 
fields. However, since the development of the Dyson-Schwinger (DS) approach to 
quark-hadron physics EE 3 the presently advocated method becomes even more 
convincing. Namely, the DS approach clearly shows how the light pseudoscalar 
mesons simultaneously appear both as quark-antiquark ( qq ) bound states and as 
Goldstone bosons of the dynamical chiral symmetry breaking (DySB) of nonper- 
turbative QCD. The solutions of Bethe-Salpeter (BS) equations for the bound-state 
vertices of pseudoscalar mesons then enter in the PVV triangle graph instead of the 
point 575 coupling, and the current algebra result © is again reproduced exactly 
and analytically, which is unique among the bound-state approaches. That the (al¬ 
most massless) pseudoscalars are (quasi-)Goldstone bosons, is also a unique feature 
among the bound-state approaches to mesons. 

A reason why the P-coupling method is simpler both technically and conceptu¬ 
ally is that the PVV triangle graph amplitude is finite , unlike the AVV one, which 
is divergent and therefore also ambiguous with respect to the momentum routing. 
Reference 1111 gives a more detailed discussion of the P-coupling method and why is 
that it is equivalent to the anomaly calculations, and illustrates this on the exam¬ 
ples of the famous decay 7r° —> 77 and processes of the type 7 —> 7r + 7r°7r~. Also, the 
PVV quark triangle amplitude leads to many (over 15 ) decay a mplitudes in agree¬ 
ment with data to within 3 % and not involving free parameters | 12 | 13 |M| This will 
be elaborated in more detail in Sec .|3 Additional advantages of this method is that 
its treatment of the 77-77' complex and resolution of the Ua(1) problem, goes well 
with the absence of axions (which were predicted to solve the strong CP problem 
but have not yet been observed and with the arguments of Ref. El that there 
is really no strong CP problem. All this will be discussed in Sec. 0 We state our 
conclusions in Sec. EJ 


2. Processes Going Through the Quark Triangle 

In this section we calculate the amplitudes for a number of processes using the quark 
triangle graphs. Figures |T| and 0 show three such PVV processes. First we consider 
7T° —> 77 decay via the u and d quark triangle graph for 7r° = (uu — dd)/\/ 2 , N c = 3 
and Goldberger-Trieman relation leading to the pion decay constant: = m/g^qq. 
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Fig. 3. Two examples of the PVV triangle graphs where just one of the vector vertices couples 
to a photon, whereas the other couples to a vector meson. These two graphs describe the decays 
of u and p mesons into a photon and a pion. 


This amplitude is finite and for th e ex perimental value of the pion decay constant, 
/„■ = (92.42 ± 0.26) MeV ^-4 givesIt2l the chiral-limit amplitude 0 of magnitude 


|^ m „=0(7T 0 -> 27)| = 


47T 2 / 7r 


= 0.0251 GeV 


very close to experimental data 


nn 


( 3 ) 


| -Texp (tT —> 27 )| — 


647 rr( 7 r° —> 77) 


m-i 


1/2 


= (0.0252 ± 0.0009) GeV 


-1 


( 4 ) 


Likewise, the u, d quark triangles for p —» 777 decay give 


□a 


I F(p 7T 7 )| = 5 %- = 0.206 GeV - 1 

/ tj - 

for g p = 4.965 ± 0.002 found from p° —> e _ e + decay^l 


( 5 ) 


r(p° 


pirn 

=-e + ) = ——| = (7.02 ± 0.11) keV . 


1277 g'l 


The calculated | F(p —> 777 )| is also near data^D 

12irT(p —> 777 ) 1 ^ 2 


F exp {p —> 777)| = 


= (0.225 ±0.011) GeV 


-1 


( 6 ) 

( 7 ) 


where q = (m 2 — m 2 )/(2m p ) is the photon momentum. [Actually, the above value 
is a weighted average of F exp (p° —> 77 ° 7 ) and F exp (p ± —> 77 ± 7 ) amplitudes.] 

Next we predict the u. d quark triangle amplitude for co —> 777 taking lo as 99% 
nonstrange^l (cos 2 </>y ~ 0.99) 


I F(u -*• 777)| 


cos 4>y e g u 


-1 


= 0.705 GeV 


( 8 ) 
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for g u = 17.06 ± 0.28 found from ui —> e e + decay. The mixing angle is : 


1 


i(4? 


K * 


- m 2 p ) - ml 


(j) V = 0 V - arctan(—p) = arctan < / - 7 —— -— 

V2 V - U 4 m K* ~ m P ) 


— arctan(— j=) 


= (5.208 ±0.092)° . (9) 

Again this theory in Eq. © is near data (0.722 ± 0 . 012 ) GeV _1 El 

Other PVV photon decays involve the p and rf mixed non-strange and ss pseu¬ 
doscalar mesons. Again the quark triangle amplitudes are a close match with data 

11 snail 4 i 

The quark -triangle (QT) calculation gives reliable predictions also for the rj and 
rf two-photon decays: 


\F(r] —> 77)1 = ' —A(5 cos (j) P — \[2 —sin 4> P ) = 0.0255GeV“ , (10) 

47H/tt 9 m s 

TTT. 

\F(p' 77)1 = TV—-p(5sin0p + v^2—cos (j> P ) = 0.0345 GeV -1 . (11) 

47H/tt 9 m s 

This should be compared with the experimental data: 


l^exp(7? -+ 77)1 = 

647rT(77 — 

77)’ 

L 

-I 

F exp (v' ->77)l = 

647rT(77 / — 
3 

♦ 77) 

777°, 



1/2 


= (0.02498 ± 0.00064) GeV 


1/2 


= (0.03133 ± 0.00055) GeV 


-1 


( 12 ) 

(13) 


where r(?y — ► 77 ) = (0.5108 ± 0.0268) keV and T(p' —> 77 ) = (4.29 ± 0.15) keV. 
The ratio of the constituent quark masses is m s /m = 2 f k / f T — 1 = 1.445 ± 0.024 


for f_± = (92.4 ± 0.3) MeV and fx = (113.0 ± 1.0) MeV 1 ^* The mixing angle is 

ITM 


(j) P 


6 p + arctan(v / 2 ) = arctan 
(42.441 ±0.019)° . 


( to )), - 2 mr K ± ml)(m^ — ml) 

( 2 m 2 K - ml - ml)(ml, - ml) 


(14) 


Next, we can calculate the p° —> 777 amplitude employing the quark-triangle 
diagram, 


I F(p° -> 777)1 = -±^_3cos <j) P = 0.456 GeV -1 . (15) 

J7T 


a We use quadratic mass formulae for mesons (See, e.g., Ref. nn and earlier). However, the input 
experimental meson masses are newest, taken from Ref. ED 
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Again, this is close to the experimental data, 


I PL_ (n° 




where q = ( m 2 p — m 2 )/(2m p ) = (194.5 ± 0.4) MeV is the photon momentum and 
r(p° —> 777) = (45.1 ± 6.0) keV. A similar situation is with the r/ —> p7 amplitude, 
for which the quark triangle calculation gives 


|F(t/ - A) I = -^3sin</»p = 0.417 GeV " 1 . (17) 

J n 

The corresponding experimental value is 


F exp (i/ -> A) I = 


47t r(?/ p° 7 ) 


ui 1 / 2 


= (0.411 ±0.017) GeV , (18) 


where q = (m 2 , — m 2 ) / (2m v ') = (164.7 ± 0.4) MeV is the photon momentum and 


T(?/ —> p °7 including non-resonant 7 r + 7 r 7 ) = (60.0 ± 5.0) keV (19) 

is the experimental decay width 1 ^. 

The q —> 7 T 7 T 7 amplitude is 


|| = | 9p ™ p °^ | = 9.80 GeV” 3 (20) 

' ' m p — s 

where s = m 2 . The rj — > 7T7T7 decay width is 

T(r? -»■ ttttj) = = 56.2 eV , (21) 

where 1), = 0.98 • 10 —5 E01. This is in a good agreement with the experimental value 


T(i 7 -> TTTry) = (60.4 ± 3.6) eV , (22) 

revealing that the vector meson dominance is the main effect, while the coupling 
through VPPP quark box loop (“contact term”) contributes little. 

It is known that u> —> 3n decay is dominated by p-meson poles. The required 
u> —> p 7 T amplitude can be estimated as 

|M vmd (u; pn)\ = | F{u -> A)l ~ 12 GeV -1 , 


(23) 
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but cannot be measured because there is no phase space for this process. The 
u —> pit amplitude is more precisely defined with QL, additionally enhanced with a 
meson loop associated with sigma exchange — 2 21 


|M( w -,p^)| Q T = |^«15GeV- 1 . (24) 

The scalar amplitude M vmd (w —> 3n) is dominated by the p meson in each of the 
three possible channels 22 ) 


|M VMD (W - 3tt)| = 2g pl rir | M(oj - ptt)| 

« 1480 GeV -3 . 

Following Thew’s phase space analysis 122 we get 


1 


L m 2 — s m 2 — t m 2 — u 


(25) 


T(w- 3 tt) = l MVM °^7 ^ m l Y * = 73 MeV 

(Z7T)° 


(26) 


wher e K , = 4.57 -10 6 is used. The predicted value is close to the experimental 
valued 


T(w -> 3 t r) = (7.6 ± 0.1) MeV . (27) 

Here we have taken to as pure NS, although it is about 99% NS, since <j>v = 
(5.208 ± 0.092)° from our Eq. ©. 

In the quark -level cr-model a quark box diagram contributes to the to —> 37 r 
decay. This box diagram can be interpreted as a contact term. It is shown that the 
conta ct c ontribution is small by itself, but can be enlarged through the interference 
effect®. 

Using cj)p = (42.441 ±0.019)° from our Eq. (1141) . we predict the tensor T —> PP 
branching ratios for 0,2(1320): 

BR( 02 -^4 = (^-) 5 2cos 2 <j>p = 2 ." 6 (data 2.96 4 0.54), 

Ci2 —> K K V rK ) 

BR( 0,2 ~^T, ) = fife 1 ) 5 2 sin 2 <j> P =0.1113 (data 0.108 4 0.025), (28) 

ct2 —> KK V J 

BR( q 2 2 ^^ ) = (^) 1 tan 2 (\>p = 0.0371 (data 0.0366 4 0.0069) , 

for center of mass momenta p^ = 535 MeV, p v ' n = 287 MeV, px = 437 MeV. The 
above data branching ratios follow from 02 (1320) recent fractions G23 BR(o 2 —> 
7 ]7 t ) = (14.5 4 1.2)%, BR(o 2 -> KK) = (4.9 4 0.8)% and BR(a 2 -> t/tt) = (5.3 4 
0.9) • 10” 3 . 
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3. Comments Related to the Gluon Anomaly 

The approach using the pseudoscalar coupling is, in our opinion, also relevant 
for the effects related to the non-Abelian, “gluon” ABJ axial anomaly. Here, we 
comment on this only briefly, and direct the reader to the original references for 
details. 

3.1. Goldstone structure and rf-iq' phenomenology 

The first point concerns the 77 - 77 ' complex and the 1 / 4 ( 1 ) problem related to it. 

In the chiral limit m n = mx = m rl8 = 0, since all members of the flavor-SU(3) 
pseudoscalar meson octet are massless in this theoretical, but very useful limit. The 
only non-vanishing ground-state pseudoscalar meson mass in this limit is the mass 
of the SU(3)-singlet pseudoscalar meson 771 . This is thanks to the non-Abelian, gluon 
ABJ axial anomaly, i.e., to the fact that the divergence of the SU(3)-singlet axial 
current 

Ki x ) = , (29) 

receives the contributions from both the finite quark masses and gluon fields G^": 

3 „2 

dnA% = 2im„ U 75 U + 2im d dj 5 d + 2i m s s~/ 5 s + e^apG^G^ . (30) 

This removes the U A { 1) symmetry and explains why only eight pseudoscalar mesons 
are light, and not nine; i.e., why there is an octet of (almost-)Goldstone bosons, but 
not a nonet. The physically observed g and rf are then the mixtures of the anoma¬ 
lously heavy 771 and (almost-) Goldstone 77 s in such a way that 77 ' is predominantly 
771 and 77 is predominantly gs- This is how the gluon anomaly can save us from the 
Ua( 1) problem in principle, and the details of how we achieve a successful descrip¬ 
tion of the 77 - 77 ' complex, are given in Refs. mmmmm Here we just sketch 
some important points. The mass matrix squared M 2 in the quark basis | uu), \dd), 
|ss) is 



~ m lu 0 0 


' 1 1 X 

M 2 = M^ a +M a = 

0 m h 0 

+ 0 

1 1 X 


L 0 0 77-4 J 


X X X 2 _ 


where M 2 A is the non-anomalous part of the matrix, since m)G = toTj = m 2 and 
m 2 s = 2m 2 K — m 2 would be the masses of the respective “non-strange” (NS) and 
“strange” (S) pseudoscalar qq mesons if there were no gluon anomaly. In the NS 
sector, in the isospin symmetry limit (which is very close to reality), the relevant 
combinations are |7r°) = \uu — dd)/\/2 as the neutral partner of the charged pions 
| 7 r ± ) in the isospin 1 triplet, and the isospin 0 combination |wu + dd)/\/2. In the 
absence of gluon anomaly, but with an s-quark mass heavier than the isosymmetric 
u and d ones, g would reduce to |NS) = | uu + dd)/V2 with the mass Tons = J 7 V, 
and g' to |S) = |ss) with the mass ms = m S s■ Both of these assignments are in 
conflict with experiment. The realistic contributions of various flavors to 77 and g' 
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and their masses (i.e., the realistic 77-77' mixing) are obtained only thanks to M\, 
the anomalous contribution to the mass matrix. In M|, the quantity fd describes 
transitions \qq) —*• | q'q') ( q,q' = u,d,s) due to the gluon anomaly and X describes 
the effects of the SU(3) flavor symmetry breaking on these transitions. In Refs. tTR 
□JI23 as the first step in solving the C/a(1 ) problem, we extract 17 $, r)± masses from 
the 77, ?/ via 

m v 8 = { m v cosOp ) 2 + ( m v ’ smd P ) 2 = (572.73 MeV) 2 , (32) 

m m = ( m v smd P ) 2 + (m,- cos dp ) 2 = (943.05 MeV) 2 , (33) 

where dp = <j>p — arctan(-\/2) = (—12.295 ± 0.019)°. The mesons r/s and r]i are 
defined as 


\Vs) =-^(\uu) + \dd) - 2\ss)) , (34) 

|t?i) =-^=(|um) + |dd) + |ss)) . (35) 

The rjs meson mass GH m Vs = 572.73 MeV is 4.56% greater than the observed ^ 
m v = (547.75 ± 0.12) MeV. The singlet 771 mass (I3H1) m m = 943.06 MeV is only 
1.56% below the observed m' v = (957.78±0.14) MeV and close to the nonstrange-ss 
mixing Ua{ 1) mass dictated by phenomenologyH SJ X S l. 2.51 


™ Ua ( l) = (3 /?) 1/2 


3 {m 2 , 

4 


- m 2 )(m 2 - m 2 ) 
m 2 K - 


= 915.31 MeV , 


(36) 


(This is also close to 912 MeV, which is the mass found in the analogous DS approach 

IraBSI ) 

We call the quantity (® the “mixing C/a ( 1 ) mass” since the mass matrix (which 
is especially clear in the nonstrange-strange quark basis) reveals that mp A (i) induces 
the mixing between the nonstrange isoscalar (|u 7 t) +1 dd) /y/2 and ss quark-antiquark 
states. Equivalently, m UA ^ can be viewed as being generated by the transitions 
among the uu, dd and ss pseudoscalar states; via quark loops, these pseudoscalar 
qq bound states can annihilate into gluons which in turn via another quark loop 
can again recombine into another pseudoscalar q'q' bound state of the same or dif¬ 
ferent flavor. The quantity j3 appearing in Eq. (|36|) is then the annihilation strength 
of such transitions, in the limit of an exact SU(3) flavor symmetry. (The realis¬ 
tic breaking of this symmetry is easily introduced and improves our description 
of the 77 - 77 ' complex considerably.) The “diamond” graph in Fig. Ogives just the 
simplest example of such an annihilation/recombination transition. Since these an¬ 
nihilations occur in the nonperturbative regime of QCD, all graphs with any even 
number of gluons instead of just those two in Fig. can be just as significant in 
annihilating and forming a C + pseudoscalar qq meson. Indeed, this nonperturba¬ 
tive Ua( 1) mass scale, Eq. 13611 . is still 3 times higher than the gluon “diamond” 
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Fig. 4. Nonperturbative QCD annihilation of a quark-antiquark bound state illustrated by the 
diagram with two-gluon exchange. The qq pseudoscalar P is coupled to a quark loop, whereby 
it can annihilate into gluons which in turn recombine into the pseudoscalar P' having the flavor 
content q'q'. 


graph evaluated perturbatively "F Thus, we cannot calculate (3 = mfj A ^/3 and 
the situation is much more complicated and less clear than in the Abelian case, as 
explained in T Can it then be founded to think that the annihilation graphs with 
the pseudoscalar meson-quark coupling, such as the “diamond” graph in Fig. 13 give 
rise to the anomalous term (3 g 2 /32'K 2 )e f _ lva pG^ v G^ in the divergence lt30l) of the 
SU(3)-singlet current A^{x), and thus ultimately to the large mass of 770 (and of 
the observed 77')? Well, this conjecture may remain a speculation since we cannot 
calculate /3 due to the nonperturbative nature of the problem. Nevertheless, when 
we use it in our approach as a parameter with the value given by Eg . 13611 . we 
obtain a very good description of the 77-77' complex phenomenology 
This includes not only the masses of 77 and 77', but also their 77 decay widths, and 
the mixing angle dp ss —13° consistently following from the masses and 77 widths. 
This gives a strong motivation for the above conjecture. 


3.2. Taming of strong CP problem 

We should also note that our c onje cture in the previous subsection goes well with 
the arguments of Banerjee et ai.*“ that there is really no strong CP problem. They 
find that one does not need vanishing © e ff = 0 — trln M. q (where M q is the quark 
mass matrix). Thus, one does not need any fine-tuning, and all CP violation in the 
QCD Lagrangian can be avoided by having 0 = 0 in its CP-violating term 

C & = -Q-£^e tlva g G ^Cf. (37) 

This term in the QCD Lagrangian breaks the J7 a( 1) symmetry and corresponds to 
the anomalous term oc /a gG^ v G^ in the divergence ll30l) of the singlet current. 
The term (pT7|) is allowed by gauge invariance and renormalizability, but apparent 
nonexistence of the strong CP violation, and also of axions, is the solid reason to 
have it vanishing. Our conjecture, that P-coupled annihilation graphs reproduce the 
effect of the gluon ABJ anomaly, naturally agrees with the vanishing of this term 
and with putting the case of the strong CP problem to rest a la Banerjee et al. . 
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4. Summary/Discussion 

We have presented and surveyed in detail the method of pseudoscalar coupling of 
pseudoscalar mesons to the “triangle” and “box” quark loops. We have reviewed how 
this method gives the equivalent results to the anomaly calculations. The P-coupling 
method has also been illustrated on the example of many decay amplitudes. 

The AVV anomaly ^ — involves 10 invariant amplitudes (reduced to 1 or 2 ampli¬ 
tudes for 7T° —> 77 decay using additional Ward identities). If instead one considers 
the PVV transition with a pseudoscalar coupling, then the P W quark triangle am¬ 
plitude is finite and leads to many decay amplitudes (ov er 1 5) then in agreement 
with data to within 3% and not involving free parameters^. To solve instead the 
former AVV decay problem, very light axion bosons have been predicted but have 
not yet been observed ^3. 

Also, there is the Ua( 1) and 0 problem involving gluons whereby strong in¬ 
teraction QCD leads to CP violation, definitely a “strong CP problem” because 
CP violation is known to occur at the 10 ” 3 weak interaction amplitude level 033. 
Physicists have tried to circumvent this u Ua{ 1) - strong CP problem” either via 
the topo logy of gauge fields or by investigating the 0-vacuum for this strong CP 
problem “l 

In this paper we have circumvented the need to deal directly with the above 
photon or gluon AVV anomalies by studying instead (finite) PW quark triangle 
graphs. Then we have given our phenomenological results - which always are in 
approximate agreement with the data. Next we return to the Ua( 1) problem and 
again use quark triangle diagrams coupled to 2 gluons. Invoking nonstrang e-strange 
particle mixing, the predicted Ua{ 1) mass is within 3% of data^- ■ ■ -I 

Thus we circumvent both photon and, admittedly on a much more speculative 
level, also the gluon ABJ anomaly without resorting either to unmeasured axions 
or to a strong CP violating term in the QCD Lagrangian. 
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